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We consider a Periodic Anderson Model (PAM) with a momentum-dependent inter-band hybridization that is 
strongly suppressed near the Fermi level. Under these conditions, we reduce the PAM to an effective low-energy 
Hamiltonian, i/eff, by expanding in the small parameter Vb/t ( Vb is the maximum inter-band hybridization 
amplitude and t is the hopping integral of the broad band). The resulting model consists of a t — J /-band 
coupled via the Kondo exchange to the electrons in the broad band, iifeff allows for studying the /-electron 
delocalization transition. The result is a doping-induced Mott transition for the /-electron delocalization, which 
we demonstrate by density-matrix renormalization group (DMRG) calculations. 
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A recent de Haas-van Alphen (dHvA) experiment on Celna 
revealed small /-character hole-pockets 1 1 1 that coexist with 
local-moment antiferromagnetism (AFM) |2|. This observa- 
tion defies the conventional view of the heavy fermion ma- 
terials according to which large-moment magnetism exists 
only for strongly localized f-electrons|3|. Therefore, a new 
paradigm for the interplay between f- and conduction elec- 
trons needs to be developed. 

The Periodic Anderson model (PAM) is the minimal Hamil- 
tonian for describing actinide and lanthanide based com- 
pounds. The model includes a periodic array of strongly in- 
teracting /-orbitals coupled to a broad conduction band via 
a hybridization amplitude Vk- Inter-band charge fluctuations 
are significantly suppressed when the bare /-electron level, 
6/, is well below the Fermi energy of the conduction band. 
In this localized regime, the PAM can be reduced to a Kondo 
lattice model (KLM) by means of a Schrieffer- Wolff trans- 
formation |4|. The /-electrons act like local magnetic mo- 
ments that can either order, or "dissolve" into the Fermi sea of 
the conduction band, leading to a large increase of the quasi- 
particles effective mass (heavy fermion phase) |5|. When e/ 
gets closer to the Fermi level, /i, the charge transfer between 
the bands increases and finally leads to /-electron delocaliza- 
tion. In general, the Schrieffer- Wolff transformation cannot 
be extended to this mixed-valence regime, when | Vk| becomes 
comparable to |e/ — Consequently, the lack of a control 
parameter poses a challenge for describing the crossover be- 
tween the localized and mixed- valence regimes. 

Here we consider a particular case that allows for extend- 
ing the Schrieffer- Wolff transformation to the mixed-valence 
regime. The basic assumption is that Vk cancels at the cross- 
ing points between the two bands. Under this assumption, 
we derive under control a low-energy effective model that ex- 
tends and generalizes the KLM. It consists of a t — J model 
in the /-band coupled via Kondo exchange to the broad-band 
electrons. The new model describes the continuous crossover 
between the localized and mixed valence regimes. We show 
that there are two general classes of the low-energy behavior 
of the PAM, depending on whether the inter-band or intra- 
/-band coherent charge fluctuations dominate. The former - 
"Kondo" - regime is characterized by strong hybridization be- 



tween the conduction and the / bands and can lead either to 
AFM or heavy fermi liquid ground state |3|. In the latter - 
"Mott" - regime, the f -electron delocalization is very similar 
to the doping-induced Mott transition. The role of the con- 
duction band is primarily limited to providing a charge reser- 
voir for the correlated /-band. Consequently, the /-electron 
delocalization is accompanied by a change in the Fermi sur- 
face topology (Lifshitz transition). Unlike the standard Kondo 
regime, new /-character pockets emerge at the transition and 
coexist with the large Fermi surface of the broad band. Ac- 
cording to calculations which are controlled only in the zero 
doping limit |[6l|3, these /-character pockets seem to be the 
characteristic Fermi surface of lightly doped Mott insulators 
(Mis). The heaviness of the /-pockets is caused by magnetic 
frustration of the kinetic energy lO . This scenario provides a 
possible explanation for the recent dHvA in Celna 1 1 1. 

We start by considering a PAM Hamiltonian of the form: 
H = Ho^ Hi, with 

Ho = -tY^ {c]^ci^ + H.c.) -l^Yl 

jo j 

Hi = ^F,,(4/z.+H.c.), (1) 

j,i 

where (j, /) indicates that j and / are nearest-neighbor sites, 
fj^ (Cj(j) creates an / (c)-electron with spin a on site j, 

^ja = fjafjo^ and rij^ = c]^Cj^. The chemical potential /i 
controls the total number of electrons. The /-electrons inter- 
act via an on-site Coulomb repulsion U. Hq contains the terms 
that do not mix the two bands, while Hj is the inter-band 
hybridization. We assume a d-dimensional hyper-cubic lat- 
tice of unit cells containing a c- (broad band) and an /-orbital 
each. The hopping, of amplitude t, is only between nearest- 
neighbor sites. This gives the dispersion = — 2t7k, with 
7k = (iCos/Cj^. Unless stated otherwise, we will as- 

sume that the inter-band hybridization amplitude is non-zero 
only between c and / orbitals that belong to nearest-neighbor 
unit cells: Vji = VoS\j.^_j..\^a is the lattice parameter). In 



momentum space we have, 

= 2yo7k , V^j = e-^"^"-^ (2) 

where is the number of /-orbitals and Vkj is the hybridiza- 
tion ampHtude for the term c^^fj^ (4^ = ^ Y.j e'^""^ c]^). 

will also assume that e j = unless stated otherwise. 
In this case, cancels at the band crossing points, and for 
U ^ 2dt, the Schrieffer- Wolff transformation 1 4 1 becomes an 
expansion in powers of the small parameter | Vk/(ek — e/) | = 
|Vb|/|t|. The resulting effective Hamiltonian is (the detailed 
derivation will be presented elsewhere | 8 1): 

kcr ka 

- '^Y.^Srs,i--/^i,i), (3) 

where Sj = \ Xlss' fjs^ss'fjs' /-electron spin on site 
J' = ^ Z)ss/(c]5^^ss'Qs' + cl^cTss'Cjs') is the conduction 
electron "bond" spin (a are the Pauli Matrices), and tji = 
T.S i^ls^is -^4s^js)' The constrained operators /j^ = /ja (1 - 
nj^) do not allow for double occupancy of the /-orbitals. f^^ 
creates an electron in the /-band with well defined momentum 

^- ^ka - ^ l^j^ Jja- 

The original broad band dispersion, ek, is renormalized to 
^k = H — ^ = ek — 2-^7k. The /-electrons acquire an 
effective dispersion e/k = e/ — 2t/7k, where = — l^o^/t 
is the effective hopping between nearest-neighbor /-orbitals 
(here we neglected terms 0{Vq /U)). The super-exchange in- 
teraction J = AP^/U is induced by the hopping If. 

Our i^eff is an extension of the standard KLM, which is 
also obtained by applying the Schrieffer- Wolf transformation 
to the PAM (H . The important difference between that well- 
known derivation and the one presented here is our original 
assumption of suppressed hybridization at the band crossing, 
which leads to a small control parameter for the perturbative 
expansion. The usual derivation of the KLM | 4 1 assumes that 
the bare e/ is below the bottom and e/ + /7 is above the top 
of the broad band to guarantee that the expansion parame- 
ter max [|Fk/(ek - e/)U^k/(e/ + /7 - ek)|] < 1. This as- 
sumption immediately implies that the /-electrons are local- 
ized and, by construction, far from the mixed- valence regime. 
We point out, however, that under certain conditions this as- 
sumption is unnecessarily restrictive. In the case we explicitly 
consider, the ratio | Vk/(ek — which controls validity of 
expansion, remains small as long as \ Vq\ <C \t\. Consequently, 
i^eff remains valid all the way from the localized to the mixed- 
valence regime. 

The qualitative picture introduced by Doniach |3 1 for the 
KLM argues that the /-moments will order antiferromag- 



netically if the Ruderman-Kittel-Kasuya-Yosida (RKKY) ex- 
change interaction between local moments is bigger than the 
Kondo temperature ksTx- For local Kondo interaction, J^, 
the RKKY exchange can be obtained by perturbatively in- 
tegrating out the itinerant electrons, JRKKY{\Yj — — 
J\ J (ike*^'^'^^ ~'^^^x(k) 13, with the itinerant static spin sus- 
ceptibility x(k) = Eq [^F(ek+q) - nF(eq)]/(ek+q - eq). 
In particular, for a half-filled band on a hypercubic lattice, the 
susceptibility diverges at the AFM wave-vector, strongly fa- 
voring AFM ordering of the local moments. In a similar way, 
the c-degrees of freedom can be integrated out in i^eff, Eq. ([3]); 
however, the result is different due to the non-trivial momen- 
tum space structure of the Kondo interaction in i^eff- It is 
easy to show, that in this case the spin susceptibility in the ex- 
pression for RKKY exchange has to be replaced by x(k) = 

Eq(7k+q +7q)^[^F(ek+q) " (Cq)] /(ek+q - Cq). Thc 7 

phase factors eliminate the divergence at the AFM wave- 
vector (in fact, xO^afm) = for half-filled conduction 
band), which makes the RKKY interaction effectively short- 
ranged and peaked at /c = 7r/2. The change in the spatial 
decay power law of the RKKY interaction (1/r^ ^/r^^^) 
results from frustration of the Kondo exchange at the Fermi 
level: the usual logarithmic divergence of x(k) at /c = 2k f 
is replaced by a logarithmic divergence of d^x{^)/dk^ at 
k — 2k F- This frustration also suppresses the short range 
amplitudes of Jrkky' Jrkky{o) ^ O.OOSFo^/t^ and 
J RKKY (2a) Fianally, the /-electron low- 

energy sector of i^eff is approximately described by a t — J 
Hamiltonian, Ht-j, in which the exchange interaction has 
two contributions: the AFM super-exchange J ^ Vq /Ut^ 
and a short range Jrkky ^ ^o/t^- Because of the small nu- 
merical prefactor in Jrkky , it remains smaller in magnitude 
than J for /7 > 2dt. 

The description of / electrons in terms of an effective t — J 
model has many important consequences. First, it implies that 
/-electron delocalization is induced by doping. For /i > 0, 
there is one /-electron localized on each orbital and the cor- 
responding moments interact via exchange. This means that 
the /-electrons behave like a MI in the strong U limit. For 
/i < 0, a fraction of /-electron density is transferred to the 
broad band. Although the kintetic energy of the /-holes com- 
petes against the magnetic ordering, the AFM correlations 
must survive for a small enough concentration of /-holes. The 
simple reason is that the kinetic energy per site scales like 
6 (concentration of /-holes) while the magnetic energy per 
site is proportional to (l-S). This leads to a phase in which 
the /-electrons are simultaneously delocalized and magnet- 
ically ordered, with an ordered moment comparable to the 
full moment. The Fermi surface of this phase includes small 
/-character hole-pockets that are disconnected from the big 
Fermi surface of the broad band. 

In the following we present numerical results computed 
with the original and the effective Hamiltonians in d = 1. 
We use the density-matrix renormalization group (DMRG) 
method 1 10] to obtain the ground state properties of both 
Hamiltonians in chains of L = 20 unit cells. These calcu- 
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lations have a double purpose. First, we show that the low- 
energy spectrum of the PAM has two qualitatively different 
regimes in the mixed-valence state. The first and most tradi- 
tional "Kondo-like" regime takes place when the average hy- 
bridization amplitude over the Fermi surface is much stronger 
than the effective hopping between /-orbitals: |Vfe^ | \if\. 
This regime is dominated by coherent inter-band charge fluc- 
tuations. In the second "Mott-like" regime, the low-energy 
physics of the PAM is controlled by coherent intra-band 
charge fluctuations, i.e., the /-electrons are well described by 
a single-band model as it can be inferred from our derivation 
of i^eff. This regime can be stabilized for \Vkp\ < \if\. The 
other purpose is to verify that i^eff provides an accurate de- 
scription of the low-energy spectrum of H as long as Vk van- 
ishes at the crossing points. 

There are several qualitative differences between the Kondo 
and Mott-like mixed- valence regimes. For d = 1, one of these 
differences appears in the momentum dependence of the /- 
magnetic structure factor S{q) — t z ' ^ 
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The Kondo-like mixed-valence regime contains short-range 
antiferromagnetic fluctuations that lead to a wide peak around 
q = TT. In contrast, in the Mott-like mixed- valence state, 
if the nearest-neighbor AFM interaction dominates (| J| > 
\Jrkky\), the /-holes carry an AFM anti-phase boundary 
ifTTI . This phenomenon is a direct consequence of the the 
intra-band nature of the /-charge fluctuations: by carrying an 
anti-phase boundary the /-holes preserve the antiferromag- 
netic alignment of nearest-neighbor moments when they hop 
between different / orbitals (see Fig[T]). Consequently, the 
peak in S{q) shifts to 7r(l ± (5), where S = 1 — Uf. 
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FIG. 1: Each hole carries an anti-phase domain wall for the anti- 
ferromagnetic correlations when a one dimensional Mott insulator is 
doped away from half-filling. 

We first analyze the Kondo regime \Vkp\ ^ \if\ in which 
the mixed-valence state is characterized by coherent inter- 
band charge fluctuations. For this purpose, we consider a 
PAM with on-site hybridization Vi^ = Vq. Figure [2] shows 
the DMRG results for different values of e/. S{q) exhibits a 
rather sharp peak at = tt in the localized regime Uf 1 
(cf = —1). The transition to the mixed- valence regime takes 
place around e/ = 0, i.e., rif becomes significantly lower 
than one for e/ > 0. The results show that the maximum S{q) 
remains at g = tt, but the peak becomes broader in the mixed- 
valence state. While the dominant magnetic fluctuations are 
still peaked at q = n, the effect of the coherent inter-band 
charge fluctuations is simply to reduce the correlation length 
of the still-dominant antiferromagnetic correlations. 

The Mott mixed-valence regime can be stabilized under 
the condition: iV^i^l <C \if\. To study this regime, we use 
the hybridization term introduced in Eq.([2]) for which we de- 



— r- 



— r- 



On-site Hybridization 
Vp/t = 0.2 
U/t = 4.0 
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FIG. 2: Magnetic structure factor S{q) for the PAM, if, with on- 
site hybridization, Vk = Vb- The results are shown for different e/ 
values, thus different different /-electron densities n/. 



rived i^eff under control. Figure [3] shows a comparison of the 
S{q) results computed with H and i^eff for different values 
of the total number of electrons Nq. Note that i^eff repro- 
duces the magnetic structure factor obtained with the PAM. 
This is indeed the expected result because the control param- 
eter Vo/t = 0.1 is small enough to guarantee the validity of 
our perturbation theory. Again, S{q) exhibits a single peak at 
the AFM wave- vector q = n in the localized state (/i > 0). 
However, the single peak splits into two symmetric peaks lo- 
cated at = 7r(l ± (5) as the system enters the mixed- valence 
regime: /i ^ and n/ < 1. As explained above, this is a clear 
signature of dominant intra-band coherent charge fluctuations. 
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FIG. 3: Comparison of S{q) obtained with the PAM and the effective 
Hamiltonian, H (dashed lines) and ifeff (full lines). Vo/t = 0.1 is 
the small parameter of the perturbation expansion that leads to i/eff • 
Different colours correspond to different values of Ne. L = 20 unit 
cells, U/t = 4.0. 



To test the relevance of the Kondo exchange term in JYeff, 
we also compare the S{q) curve obtained with the PAM (for 
Vk given by Eq.(|2])) against the results for a pure t — J model, 
which neglects the Kondo coupling present in Eq. ([3]). The 
comparison is shown in Fig.Qa) for Fo/t = 0.2. The good 
agreement confirms that t— J, as well as the underlying single- 
band Hubbard model, provides an accurate description of the 
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/-electrons in the mixed-valence regime under consideration. 
As we discuss below, this fact has important implications for 
more realistic higher-dimensional systems. To test the robust- 
ness of the double-peak structure in S{q) away from the per- 
turbative regime, we also compute S{q) for the PAM with Vk 
given by Eq.([2]) and Vo/t = 0.5. The results are shown in 
Fig. |4|b). Although the peaks become broader, the double- 
peak structure remains robust in the mixed- valence regime. 
This indicates that the mixed state is still dominated by coher- 
ent intra-band fluctuations, even away from the perturbative 
regime \Vo\ <C \t\. 
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FIG. 4: (a) Comparison of S{q) obtained with the PAM, H, (full 
lines) and a pure t — J model for the /-electrons (dashed lines) for 
Vo/t = 0.2 and e/ = 0. (b) S{q) for the PAM with Vo/t = 0.5 and 
e/ = 0. The double-peak structure remains robust away from the 
perturbative regime, (c) Same comparison as in (a), but for e/ = 0. It 
and the same Uf for both models. L = 20 units cells, U/t = 4.0. 

Here we considered a particular form of Vk [Eq.(|2])] in order 
to prove that the mixed valence regime can be dominated by 
single-band physics, even starting from a completely flat bare 
/-band. In general, the hybridization will not cancel exactly 
at the Fermi level. As an example, in Fig. |4jc) we show a 
comparison between the t — J and PAM models for similar 
values of n/ and e/ = O.lt. This value of e/ leads to a non- 
zero hybridization at the Fermi level comparable to the / — / 
effective hopping: \Vkp\ — \if\/2. Although the double peak 
structure of S{q) gets broadened forn/ = 0.9, it remains well 
defined forn/ = 0.8 and 0.7; this indicates the robustness of 
the Mott mixed- valence regime proposed here. Furthermore, 
a realistic PAM should also include a bare / — / hopping t / . 
The single band physics derived in this work remains robust 
as long as the average hybridization over the Fermi surface 
remains smaller than the effective / — / hopping: |Vfc^| <C 
\if -\- tf \ |8|. This observation extends the relevance of our 
results beyond the particular PAM considered in this work. 



The physics of lightly doped Mott insulators is well under- 
stood in the zero concentration limit (one hole) |7 1. As long as 
the system remains antiferromagnetically ordered, the quasi- 
particle bandwidth is of order of the effective exchange be- 
tween local moments. Each quasi-particle consists of a hole 
(/-hole in our case) dressed by a local antiferromagnetic dis- 
tortion. The effective mass of the magnetic distortion can be 
much bigger than the mass of the bare hole, in which case 
the effective quasi-particle mass m* is dominated by the ex- 
change interaction, m* ex 1/J (in d > 1) ITJ. In this way, 
heavy fermion behavior can originate from and coexist with 
local-moment antiferromagnetism (ordered moment compa- 
rable to the full moment) in systems with no more than one 
/-electron per ion, such as the Ce-based compounds. 

Our work provides a scenario for the f-electron delocal- 
ization that explains several qualitative aspects of a recent 
dHvA experiment in Celns (H, which cannot be accounted 
for within the conventional view of the heavy fermion materi- 
als. In addition, coexistence of local moment AFM and heavy 
electron superconductivity was observed in the related layered 
compound CeRhlns under pressure |[l2l[T3l. This is another 
strong indication that the heavy-fermion behavior can coexist 
with the local moment AFM. Future experiments are expected 
to clarify the applicability of our scenario to these and other 
heavy fermion compounds. 
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